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ABSTRACT 

A Lie-algebraic approach successfully used to describe one-dimensional Bloch oscillations in a 
tight-binding approximation is extended to two dimensions. This extension has the same alge- 
braic structure as the one-dimensional case while the dynamics shows a much richer behaviour. 
The Bloch oscillations are discussed using analytical expressions for expectation values and 
widths of the operators of the algebra. It is shown under which conditions the oscillations sur- 
vive in two dimensions and the centre of mass of a wave packet shows a Lissajous like motion. 
In contrast to the one-dimensional case, a wavepacket shows systematic dispersion that depends 
on the direction of the field and the dispersion relation of the field free system. 
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1 Introduction 

The quantum dynamics of a particle in a tilted periodic potential, a so called Wannier-Stark 
system, has been the subject of investigation from the beginning of quantum mechanics [1-3]. 
In the early years the motivation was to understand the properties of an electric current in solid 
crystals while in recent years the realisation of semiconductor superlattices and the experimen- 
tal progress in manipulating cold atoms resp. Bose-Einstein condensates in optical lattices has 
motivated new theoretical studies in this field (see [4] and the references given there). Bloch 
oscillations are one of the striking phenomena in such systems and an example for a counter- 
intuitive behaviour in quantum mechanics [4-11]. While [4-7] discuss the one-dimensional case 
analytically in a single-band tight-binding model, in [8-11] two-dimensional Bloch oscillations 
are considered, however, only partially in an analytical treatment. 



The aim of the present paper is to supplement the discussion of Bloch oscillations in two- 
dimensional lattices by an analytic treatment within the tight-binding model. The Lie-algebraic 
approach to the tight-binding model in [5] is extended to two space dimensions using the tech- 
niques of algebraic structure theory. Within this model, not only the nearest neighbour couplings 
along the axes (the so called von-Neumann neighbourhood) is taken into account, but also the 
interactions along the diagonals (Moore neighbourhood) will be considered and explicit expres- 
sions for the expectation values of the position operators will be derived. It is shown that Bloch 
oscillations in two dimensions occur only for special initial conditions and that typically the 
oscillation is superposed by a directed motion orthogonal to the field. On the basis of this ex- 
amination the results can be extended to the single-band model where all couplings within one 
band are taken into account. 

A two-dimensional Wannier-Stark system can be described by a Hamiltonian with a periodic 
potential V and an additional linear potential: 



and r n = n\e.\ + n^e-i, Ui £ Z. For simplicity the period in both directions is set equal to one. 
The force F = Fo+F u (t) consists of a static part and a part with an arbitrary time dependence. 
A convenient approximation is to expand the Hamiltonian in a basis of Wannier-functions \n, I) 
which are localised at site n = 712) of the periodic potential. The index I labels the different 
bands. Taking into account only a single band and the coupling between adjacent sites, the 
Hamiltonian for the one-dimensional system can be written as 



where \m) are the one-dimensional Wannier-states and A is the band width of the energy 
dispersion relation for the field-free system. An alternative tight-binding approach using an 
expansion in the eigenstates of the Wannier-Stark Hamiltonian (^Q) for a constant field can be 
found in [12]. 

The tight-binding Hamiltonian can be easily extended to two dimensions using the two- 
dimensional Wannier-states \m,n) with (m, n\m', n') = 5 m * m S n / n . Then, the Hamiltonian reads 
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Introducing the operators 
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and 

+00 +00 



K\ = \m, n) {m + 1, n\ , K 2 = \m,n){m, n+ 1| (5) 



m,n=— 00 m,n=— 00 



together with their adjoint operators k\ and if] the notation can be simplified: 

H = ^{K 1 + K\) + ^-(K 2 + k\) + F l N l + F 2 N 2 . (6) 

The operators K\ and i^2 are unitary, i.e. they fulfil the relation 

KT x = k) for j = 1,2 (7) 

and describe the coupling of adjacent sites along the two axes. The hermitian operators N± and 
N 2 can be interpreted as 'position' operators in direction 1 resp. 2. 

For later convenience, the variable / = F/h will be used instead of the field F = (Fi,F 2 ) 
and it is assumed that the components f\ and f 2 of the field / have a rational ratio, i. e. 

= ~ with q, r € Z coprime . 
h r 

Then, the field can be written as 

fh\ f M 



J2) Vq 2 + 



r 



(8) 



with / = |/|. We will assume that the direction of the field is constant in time, whereas its 
magnitude / can be time-dependent. For the example of Bloch oscillations studied below, the 
field is constant in time. 



2 The Lie-algebra 

The commutators of the operators Nj and Kj can be directly calculated using equations (0} and 
(jSJ) and the orthogonality of the Wannier-functions. The non-vanishing commutators are 

[k j ,N j ]=K j , [K],N j ] = -k] fori = 1,2. (9) 
So the operators {N±, K\, k\, N 2 ,K 2 , K^} form a Lie-algebra 

c = d e £ 2 = {Nuk^kl} e {iv 2 , k 2 , k\) (10) 

with [£i,£ 2 ] = {0}. Since this algebra can be decoupled into two disjoint algebras which 
correspond to the Lie-algebra of the one-dimensional tight-binding system [5], this expansion 
is trivial. One can achieve further insight into the dynamics of the system by extending the 
algebra by the operators 

+00 

kxk 2 = ^2 \m,n)(m + l,n + l\, (11) 

m,n=— 00 
+00 

Kik\= ^ \m,n)(m + l,n-l\ (12) 

m,n=— 00 
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together with the adjoints k\k\ and k\K2. The ordering of these operators is arbitrary since 
they commute. Adding these operators to the system takes couplings between adjacent sites 
along the main diagonals into account and, from the algebraic point of view, couples both 
algebras C x and C2 ■ The non-vanishing commutators for these operators are 

[k x k 2 , n x ] = k x k 2 [k x k 2 ,N 2 ] = k x k 2 

\k\k\M =-k\k\ [k\klN 2 ] =-k\k\ (13) 

[k x kl,$r x ] = k x k\ \k x k\,N 2 \ =-k x k\ 

[k\k 2 ,N x ] =-k\k 2 [k{k 2 ,N 2 ] = k\k 2 . 

Thus, this set of operators closes under commutation and forms a Lie-algebra 

c = {n x , k x , k\, n 2 , k 2 , kl k x k 2 , k\k\, k x k\, k\k 2 } (u) 

with the corresponding Hamiltonian 

H = ^(k x + kl) + ^(k 2 + kl) 

+ (k x k 2 + k\k\ + k x k\ + k\k 2 ) + nf x N x + hf 2 N 2 . (15) 

The factors k = 1, 2, 3 are the band widths of the dispersion relation in the direction of the 
axes (Ai, A2) resp. the main diagonals (A3). A discussion of the spectral properties of such 
two-dimensional tight-binding Hamiltonians taking into account various couplings can be found 
in [13]. 

For later convenience we introduce the notation 

k UjV = k%kl with u,v(£Z, (16) 

where the commutators read 

[K U)V , Nj] = e u , v j K UyV for J = 1,2; u, v G Z (17) 
with e U)V) j = u5j tX + v 5j t2 ■ Then the Lie-algebra can be written as (cf. equation ((7J)) 



£ = {N X ,N 2 ,K U , V 

with 



(u,i;)eMH (18) 



M = {-1,0,1} 2 \{(0,0)}, (19) 

i. e. the identity operator I = -Ko,o is not included in the algebra. In section [5] the algebra will 
be extended by taking into account all possible operator products. Then, the Lie-algebra has 
still the form (Q2J), but with M = I? \ {(0,0)}. 

In this algebraic approach to the tight-binding model, all calculations can be done with- 
out using the explicit representations of the operators in (J3J and ©, only making use of the 
commutation relations in equation 1)17(1 . 



4 



3 The time evolution operator 



As described in [5, 14], the Lie-algebra can be decomposed into a semidirect sum, C = 1 (ESS, 
with 

1= [k u>v (u,v) e m} and (20) 
SS = {N 1: N 2 }. (21) 

Here, X is an ideal of the Lie-algebra and SS is a subalgebra. The Hamiltonian can be written 
as a sum of operators, H = Hs + Hi, with Hs € SS and Hj S X: 

H s = ft/iiVi + hf 2 N 2 , (22) 
Hj = ^-(Ki,o + K-i, ) + ^(Ko,i + K ,-i) 

+ -j-(Ki,i + K-i,-i + Ki,-i + K-i,i) ■ (23) 
In order to find the time evolution operator U with 

iH^ = HU, (24) 
dt y 1 

one can make use of the algebraic structure, which allows a factorisation of the evolution operator 
as U = Us Uj with 

ih ^ = H S Us , m ^ = {Us'HjUs) Ui (25) 

(see [14-16] for the details). Thus, the problem of calculating the time evolution operator is 
splitted into two parts. First one has to solve the differential equation for the subalgebra SS 
and afterwards solve the problem for I. 

Before we begin with the derivation of the evolution operator, we introduce the so-called 
r-evolved operator A, defined by 

e 2adf i = e 2f ie- 2f 

= A + z[t,A]+^[t,[f,A)] + Z ^[t, [f,[f,A]]] + ... . (26) 

with f , A € C. 

For the Lie-algebra C these expression can be calculated directly using equation (|17|). The 
non-trivial relations are 

e 2 adJ ^ Nj = Nj + ze u , vJ K u , v and e 2 ad ^' K u ,v = e~ z ^k u , v . (27) 
After these preparations the first equation in (|25jl can be solved and yields 

U s (t) = e-^ 1 ^?^) (28) 
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with 



(i) 



A(r)dr 



f f(r)d 
Jo 



r = qr)t, 



m 



(2) 



/ 2 (r)dr = D]i, 



/(r)dr. 



(29) 



s/q 2 + r* Jo 



To solve the second differential equation in (|25|) . one has to evaluate the right-hand side, 



Ug^HiUs 



A 



Ah 

+ — 

4/i 



Ah 



e~ irr,t K 01 + e +irr *K { 



+ e +i( ' ? - r), ' t i^_i i i 
which is an element of X. The solution of equation (|25[) then reads 

= exp {-i[xf ' 0) ^0 + X^'^-^ + xf' l) K^ + Xf 



0.1 



(30) 



with 



and 



Xt 



M =g (u,v)f e -i(uq+vr) V r dr > ( n , w ) € 



(31) 

(32) 



o (±l,0) _ ^1 (0,±1) _ ^2 (±1,±1) _ (±1, T 1) _ ^3 

5 4/i ' 5 4ft ' 9 9 Ah' 

The function rjt is defined in equation (|29|) . Note that the function x| is the complex 



conjugate of x\ 



(u,v) _ 



(—u,—v) (u,v)* 

Xt -Xt 



(33) 



Time-independent fields: 

So far, the expressions are valid for an arbitrary, time dependent field /. For the time- 
independent case, / = const, the integration of (j2H) yields 



(i) 



,J t and rj^ 



rf 



yV + r 2 



i. 



Defining the Bloch time Tg for the two-dimensional system via 

2vr 
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(34) 



(35) 



and the Bloch frequency by lob = 2tt/Tb, the evolution operator Us(t) can be written as 



The functions Xt can a l so be evaluated directly, 



t 





2 sin ((«« + for «ff + w ^ 




(37) 

for tig + = 

Since both operators Nj have an integer valued spectrum, the evolution operator over one Bloch 
period simplifies to 

U(T B ) = Us(T B )Ui(T B ) = Ui(T B ) • (38) 

Only if the field is not directed along one axis or a main diagonal, the operator Uj(Tb) as 
defined in (|3Tj) is equal to the identity operator, Ui{Tb) = I, since then uq + vr ^ holds for 
all (it, v) € M and therefore all functions vamsn according to equation (|37|). Otherwise, 

there exists one function Xt'^ W1 ^h uq + rv = 0, namely Xt' ^ : that does not vanish according 
to xfe~ 9) = <^-«)T B . 

In other words, if the coupling in the direction orthogonal to the field vanishes, g( r ~ q "> = 0, 
resp. is not included into the model the evolution operator over one Bloch period is equal to 
the identity operator with the consequence that the dynamics of every initial state is periodic 
in time. 

In addition, it should be pointed out that in the time-independent case the evolution operator 
can be directly written as 

U(t) = exp{-iHt/h} 

= exp{ -i [g^ (K lj0 + jLi )0 ) + <7 (1,0) (#o,i + #o,-i) 

+g( 1 ' 1 )(K lil + £_! _! + K X -x + ^_l,l) + f t Nl + f 2 N 2 ] t} 

= exp{-i[ Yl g^K^v + quB^+riVB^t} (39) 

(u,v)eM 

with g(~ u ~ v ) = g( u ' v ) . The connection between this equation and the product form U(t) = 
Us(t)Ui(t) is given by the application of Taylor's theorem introduced by Sack [17]: 

exp{ K (i + AS)} = e nA exp {{\B /a){l - e~ Ka )} (40) 

where A is a shift-operator, [A, B] = aB, and A, k, a (a ^ 0) are complex numbers. This 
relation can be directly extended to an arbitrary number of commuting shift-operators B^ with 
[A, B„] = a^B^ (a„ + 0) and \B^ B v \ = 0: 

exp{ K (i + £ A^)} = e Ki n ex PU A ^K)(l " e"^)} • (41) 

Identifying 

i = giVx + WV 2 , B„ = K^, K = -iu B t, \a = g {u ' v) /u B , a^ = -qu-rv, 
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the evolution operator can be written in the product form 

'} n 



U(t) = exp<j -v J j B (qN 1 + rN 2 )tj J_ | (\xp<j <j 

(u,v 



-i(qu+rv)ui B t _ 1 
(u,v) ~ * 

(qu + rv)uj B 



K, 



u,vj ) (^) 



which is of course identical to U(t) = Us(t)Uj(t) given in equation Q28JI and (|3*Tj) , 



4 Expectation values 



In this section the time evolution of the expectation values is discussed in order to get a better 
understanding of the dynamics. The time-dependence of the expectation values is derived in 
the Heisenberg picture. Using the expressions in equation (|27jl the following calculations can be 
carried out without difficulty, however, the expressions appear a bit extensive due to the number 
of involved operators. The Heisenberg operators for N\(t) resp. N 2 {t) read 

Aq(i) = U^^NtUit) = U-^t^U^t) = 
Ni + i 



JWfr J-^K -lJ 1 ' 1 )^ J- l ~ l ^K 

Xt K i,o~Xt K -i,o + Xt K i,i~Xt 
+ Xt K i,-i~Xt 
N 2 (t) = UfWNzUi® = N 2 + i[xJ 0,1) ^o,i - xJ°' _1) ^o,-i 

+ Xt K i,i~Xt &-i-i~Xt fti-i+Xt K -i,l 



(43) 



(44) 



Operators in the Heisenberg picture are indicated in the notation by the explicit time argument. 
The only differences between the operators N\(t) and N 2 (t) are the changed sign for the term 
with the operators K\-\ and K~i t \ and the expressions xl '^ -^o,±l instead of Xt^ -^±1,0- 
More involved is the calculation of N±(t) = [Ni(t)] 2 which can be obtained by squaring equa- 
tion (1331) : 



iVf(t) = iV 1 2 + 2|xr' 0) | 2 + 2|xJ 1 ' 1) | 2 + 2|xJ 1 '- 1) | 2 



, (1,0)2 9 (1,1) (1-1) .A. 

[Xt + 2 Xt Xt ) K 2,o 



, (-1,0)2 (-1-1) (-1,1), A v (14)2^ , (-1,-1)2 A 

+ {Xt + l Xt Xt ) K ~2,0 - Xt K 2,2 + Xt X--2-2 



(1,-1)2 A (-1,1)2 A- . . 

Xt K 2,-2 + Xt + 1 



(1,0) f (-1,0) ft , v (l,l) J T> 

Xt J i ~Xt J i+ Xt Ji K o,i 



(-1,-1) flfr (1,-1) f (-1,1) 1 » (1,0) (1,1) f> 

~Xt 'h K o,-i + Xt Ji K o,-i-Xt Ji K o,i - *Xt Xt K 2,i 

(-1,0) (-1,-1) r> , „/ (1,0) (-1,-1) (-1,0) (1,-1) s A 

~ l Xt Xt K~2,-l+*\Xt Xt +Xt Xt ) K 0,-i 

I „/ (-1,0) (1,1) (1,0) (— 1,1)\ -f> 9 (1,0) (1,-1) A 

+ 4Xt Xt +Xt Xt )^o,i-^Xt Xt A 2-i 

9 (-1,0) (-1,1) a „ (1,1) (-1,1) a 9 (-1,-1) (1,-1) A 

- l Xt Xt K -2,l + *Xt Xt K 0,2 + *Xt Xt K 0,-2 • 



Here, the operator Jj is the anti-commutator of N~ and Kj, 

J J = [N J ,k J ] + = N J k,+K 3 N J 



(45) 



(46) 



The expression for N 2 (t), obtained by squaring N 2 (t), has the same structure as N'f(t). Instead 
of the functions Xt' 
K u v are exchanged. 



of the functions xi ±1 '°' > the functions Xt'^ appear and, analogously, the indices of the operators 



8 



Time-independent case: 

Now we turn back to the time-independent case in order to analyse the Bloch oscillations in 
more detail. The expectation values for the operators K u>v , Ji-ftTf and Kf J 2 can be decomposed 
as 

(K u>v ) = K UtV = \K u>v \e iK ^ (47) 
{J l K v 2 ) = J v = \J v \e^ (48) 
(K?j 2 ) = L u = \L u \j^ , (49) 

where the operator Kf resp. K 2 corresponds to K u ,o resp. Ko,v in the alternative notation Q16JI . 
The expectation value {KfJ 2 ) only appears in the expression for {N 2 )t (cf. equation (|59|l). All 
expectation values including only the unitary operators K u ,v have modulus smaller than one 
except for the eigenstates of these operators where they have exactly modulus one. 

In order to calculate the expectation values for the operators Ni(t), N 2 (t) and Ni(t) the 
functions Xt are a ^ so decomposed into their moduli and phases, 

Then, the desired expressions for the expectation values read 

{Ni) t = (iVOo + 2IXJ 1 ' )! \K h0 \ sin (<?^ - k 1i0 ) 
+ 2| x ; i ' 1) ||iv 1 , 1 |sin(( (? + r)^-K lil ) 

+ 2| x ; i '- 1) ||iv 1 ,_ 1 |sin(( (7 -r)^- Kl ,_ 1 ) , (51) 
(JV 2 ) t = (iV 2 ) + 2| X | 0,1) I \K 0A I sin (r^f* - « ,i) 
+ 2| X l 1 ' 1) ||iv 1 , 1 |sin(( (? + r)^-K 1 , 1 ) 

- 2[ X , C1 '- 1) | |*T W | sin (^ - r)^f* - (52) 



and 



(JV?) t = (JV?> + 2| X i 1 ' 0) | 2 (1 - \K 2fi \ cos (2^ - K2i0 )) 
+ 2| X f' 1) | 2 (1 - \K 2 , 2 \ cos {2(q + r)*f - n 2 , 2 )) 

+ 2\ X ^ 1] \ 2 (1 " 1^2,-21 COS (2(q - r)^f - K2,_ 2 )) 

+ 2| X f' 0) I I Jo| sin - Mo ) + 2| X J 1 ' 1) 1 1 Jx| sin ((, + r)^ - Ml ) 

+ 2| x ; i '- 1) ||J_ 1 |sin(( (? -r)^- / i_ 1 ) 

- 4| x ; i ' 0) xS 1 ' 1) | \K 2tl \ cos ((2g + r)^ - « 2il ) 
+ 4|xf'° ) xi 1 ' 1) |^o,i|cos (r^-Ko,i) 
-4| X ; i>1) xS 1 *~ 1) ||^2,o|cos(2^- K2 ,o) 

+ 4|xJ 1 ' 0) xl 1 '~ 1) ||^o,i|cos (r^-«o,i) 

- 4|x; i '° ) xS 1 '" 1) | |#2,-i| cos ((2g - r)*f - K2,-i) 

+ 4|xf' 1) xf'~ 1) | l^o, 2 |cos (2r^ - ^,2) • (53) 
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It should be pointed out again that the functions Xt" depend only on system parameters 
(g^ and /) while \K U>V \ and K UtV resp. \J U \ and fi u are determined by the initial state. As 
long as the field is not directed along the axes or the main diagonals, i. e. 5 / / r, and 
q + r^O^q — r, the expectation values {Ni)t, (-^2)*, (A^i)t and (iV|)t and therefore the 
widths A^-.(t) = {Nj)t — j = 1,2, are periodic in time since all functions Xt^ are 

periodic according to equation 1)37(1 . Thus, the results from the one-dimensional case concerning 
Bloch oscillations and breathing modes remain valid: A wavepacket performs a two-dimensional 
Lissajous-like oscillation which shows no systematic dispersion, i.e. Afj (t) is bounded in time. 

More interesting is the case that the field is directed along one axis or a diagonal, i. e. 

(q, r) = (1, 0), (0, 1), (±1, ±1), (Tl, ±1) • (54) 

According to equation (|37|) . the function xt' ^ grows linearly in time, Xt' ^ = g^'" 1 ^ t. Since 
the argument (uq + vr)ujBt/2 of the related sin-function in the expressions for the expectation 
values {Nj)t (cf. equations (|51|). I|52|l) vanishes for all times t, only the phase factor n r _ q re- 
mains. This has the consequence that the term containing Xt ^ m ^ ne equation for (N\)t resp. 
{N2)t increases linearly in time as far as the corresponding phase factor k U)V does not vanish 
either. 

Thus ({Ni)t, {N2)t) performs an oscillatory motion (due to the remaining terms) with a super- 
imposed linear motion orthogonal to the field direction. 

For (q,r) = (1,0) it is clear that the linear motion is orthogonal to the field, since the linearly 
growing function |xt°'^ I = |xi°' ^ I appears only in the expression for (iVj)*. The same argument 
holds for (q,r) = (0, 1). If the field is directed along one of the main diagonals the motion is 
also orthogonal to the field as can be checked directly. In general, the velocity of the directed 
motion is determined via equation (|27[) by the commutators (|17|) and can be written as 



« = 2^)|X- Pi _ ff | B in(« r ,_ 9 ) p- 9 ' 1 ] =2g( r >-4\K rt -. q \ B m(Kr,- q ) ( 



(55) 



Independent of the initial state, {Nf)t grows quadratically in time due to the second up to 
fourth term in equation (|53|) . The contribution from the other terms is at most linear. The case 
that the factors cancel each other can be excluded because this requires |.K" Ui „| = 1 which is only 
fulfilled if the initial state of the system is an eigenstate of the shift operator K uv . But these 
states are delocalised in a basis of Wannier-states and can therefore be excluded. With {N?) 
also the variance A^.(i) = (Nf)t — J = 1,2, grows quadratically in time. Only for states 

with special initial phases k uv and amplitudes a localisation effect can be achieved. But 

this is no systematic effect, i.e. not depending on the field and the coupling constants, but only 
caused by a specially prepared initial state. In section El this behaviour is analysed in more 
detail for a Gaussian initial state. 



Separable Case: 

Before we conclude this section, a brief discussion of the separable case, where only couplings 
along the axes are taken into account, allows further insight into the dynamics. Here, we 
also assume that the field is time-independent. For the separable system with A3 = in the 
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Hamiltonian (|15|). the expressions for the expectation values, equations (|5T]) and (|52*|) . simplify 
drastically: 

(jfyt = (iV 1 ) + 2|x| 1 ' 0) ||Ki,o|sin(^- Kl ,o) 



=° (#i>o + — [*i,o| sin (<^) sin (g^ - « 1)0 ) 
= (Ni) + 7^—\ K i,o\ [cos(qu;Bt - «i,o) - cos(« 10 )] 

2<?o; J B 

(JV 2 ) t = (iY 2 ) + 2|x| ' 1) ||Ko,i|sm(r^- Ko ,i) 



|* ,i| [cos (rw B t - « ,l) - cos(k Q) i)] 



(56) 



(57) 



If the field is not directed along one of the axes, the result is a Lissajous-like oscillation of the 
wavepacket depending on the field direction (q, r) and both initial phases k^o and Ko,i (and not 
only the phase difference Acf) = Ki t o — Ko,l)- According to equation (|52[). the expectation values 
for Nf(t) read 

(N?) t = (JV?) + 2| x r' 0) | 2 [1 - \K 2t0 \ cos {2q^L _ „ 20 )] 



+ 2| X J 1 ' 0) ||J |sin (g^-Mo) 



(58) 



(JVf ) t =(iV 2 2 } + 2|^ U ' 1) | 2 [1 - |* , 2 | cos (2 



,<£Bt 



«o,a)] 



+ 2| X J°' 1) ||Lo|sin (r 



(59) 



These are periodic functions of time if the field direction is not along one of the axes. The 
variances A 2 ^ (t) can be calculated directly: 

,(1,0) .2 



+ 2\ X r } | 2 1 - 1*2,01 cos (2q^f - K2fi ) - 2\K lfi \ 2 sin 2 (q*f - Kl , )] 



(1,0), 



I J 1 sin (q 



Ho) - |*i i0 |(iVi)osin (. 



Mi,o) 



(60) 



10,1) ,2 



+ 2|x^'T [l-|*o,2| cos (2r^f* -«o, 2 ) - 2|* 0jl | 2 sin 2 (r^ - k ,i)] 



+ 2|Xt 



(0,1), 



|L |sin 



^0 



|* ,i|(iV 2 }osin (rUf-no) 



(61) 



For a field direction along one axis, e. g. the first, a wavepacket shows dispersion in the direction 
of the second axis due to Xt = g^'^t. Besides, (N%}t can increase also linearly in time, if the 
initial state provides a phase ko,i n7T i n ^7L. 

Concluding, for a field directed along an axis or the main diagonal the Bloch oscillation of the 
centre of mass of a wavepacket is superposed by a directed motion orthogonal to the field. Only 
for special initial conditions this directed motion is not present and one finds strictly periodic 
oscillations. Nevertheless, the projection of the motion onto the field direction is always periodic. 
For other field directions the Bloch oscillations still remain, since the couplings orthogonal to the 
field, g( r ~ q \ are not included into this tight-binding model. But, extending the tight-binding 
approximation by taking into account the coupling to more directions leads to a superposed 
directed motion also in these cases, as will be seen in the next section. 
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5 Single-band models 



As in the one-dimensional case (cf. §9 in [5]) one can add an arbitrary number of additional 
shift operators 

k x k x k 2 , k x k 2 k 2 , k x k x k x k 2 , k x k x k\, ... (62) 

to the algebra. Physically these operators couple adjacent sites in the directions given by their 
indices {K u ,v with u and v coprime, is the coupling in direction (u,v)) or describe the coupling 
to further sites (K nUt1lv with n € N and u, v coprime is the coupling to the nth site in direction 
(u, v)). Generalising expression (|18|) . the Lie-algebra can be written as 



C=\N X ,N 2 ,K U , 



[u,v) 



i} with M = Z 2 \ {(0,0)}. (63) 



The corresponding Hamiltonian is a direct extension of the tight-binding Hamiltonian and reads 

H= 9 (u ' v) ku,v + fiN 1 + f 2 N 2 with g (u,v) = g (-u,-v ) (64) 

(u,v)eM 

The evolution operator can be calculated in the same way as in section |31 U(t) = Us(t)Ui(t) 
with Us(t) defined in equation (|28|) and 

= expj-i £ X^kuA. (65) 

Instead of only eight functions xi™'^ i n ^ ne tight-binding approximation, one gets infinitely 
many functions. Therefore, for every rational field direction (q, r) there exist pairs (nr, — nq), 
n £ Z, so that = g^ u,v ^ t with the same consequences as in the tight-binding approximation. 

Especially, the velocity of the directed motion is given by equation (|55|), But, one has to keep in 
mind that for realistic potentials, as investigated for example in [9], the factors g( u ' v > decrease 
exponentially with increasing |u| + so that the directed motion and the dispersion are very 
small compared to the oscillation. 



6 Case study of a Gaussian initial distribution in a Wannier 
basis 

In this section, the results from section 0] will be applied to an explicit representation of the 
initial wavepacket in a basis of Wannier functions, 

IV>(<)> =X) c "v>(*)K n >- ( 66 ) 

m,n 

The expectation values can be calculated using expressions (foTj) . (|52"|) and The expec- 

tation values (K U)V ) = K u ^ v can be evaluated as 

(K UiV ) = (i;\k?k v 2 \4>) = J2 c^ n c m+u , n+v (67) 

m,n 
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with Cm^n = c mtn (0) (cf. equation ©). In the following, we will assume that the initial distri- 
bution is Gaussian, i. e. 

c m ,n = at exp <^ 1- \(kim + k 2 n) \ (68) 

with normalisation constant a € R. Then, the expectation value ()67|) can be written as 

(K uv ) = a 2 e K k i u + k 2v) ^ e -(m 2 +n 2 + {m+u) 2 +(n+v) 2 )/(4,(T 2 ) 

= \K U:V \e i ^ u+k2V ^ (69) 

and the phase factor in equation (|1T|) is given by 

k u ,v = k\u + k 2 v . (70) 

Reconsidering the expressions for {Ni)t and (N 2 )t in equations (|51() and ()52l) . the behaviour of 
the wavepacket depends crucially on this phase K u ,v The function xf' 9 \ which describes the 
coupling strength in the direction orthogonal to the field, grows linearly in time according to 
equation (|3Tj) . Then, depending on the argument of the related sin- function, sin(K rj _ g ), this 
linearly growing term contributes to the expectation value or is cancelled by sin(fi; rj _g) = 0. If 
K r _ q is an integer multiple of tt, this term does not contribute and the expectation value is 
periodic in time. That means that for a given field direction (q, r) one has to make sure that 
the condition 

K rj -q = k\r — k 2 q = nir with ngZ (71) 

is fulfilled to get a periodic centre of mass motion. If this condition is not fulfilled, the periodic 
motion is superposed by a directed motion orthogonal to the field. The left-hand side of figure E 
illustrates this behaviour. Displayed are the expectation values {Ni) and (N 2 ) as a function of 
time which are obtained by a direct numerical integration of the time-dependent Schrodinger 
equation for the Hamiltonian (|15|) and parameter values h = 2.81/2-7T, F = 0.1, Ai = A2 = 1/4, 
A 3 = 0.008. 

Independent of the phase (|70|) of the initial wavepacket, the width of the wavepacket increases 
quadratically in time. Considering for example equation ()53j) and a field along the second axis 
or one of the main diagonals this can be seen in the following way. The relevant term for the 
quadratical dispersion has the form 

2\xi r '~ q) \ 2 (1 " \K 2r - 2g \ COS («2r,-2,)) • (72) 
Together with the relevant term for {Ni)t, 

2| X ! r '- 9) ||K r ,_ 9 |sin( Kr ,_ (? ) (73) 
one gets for the width of the wavepacket 



A^ i (t) = (iV 1 2 ) t -(iV 1 ) t 2 

(r ,-?) |2 n _ I ir„ „ l™ D ^_ „^_Ol^ |2, 



« 2\ X r T [1 - \K 2r „ 2q \ cos (K2r,_2 9 ) - 2|K r ,_ 9 | 2 sin 2 (k^,)] . (74) 

The phase factor K 2r _ 2q can be written as 2K r - q using equation (|70|) so that the following 
inequality holds 

1 - \K 2r - 2q \ COS {n 2r -2q) - 2\K r - q \ 2 SIB. 2 (/C r ,-g) 

= 1 - (2\K r _ 9 | 2 - \K 2r -2q\) sin 2 (n r _ 5 ) - |if 2 r ,-2 9 | cos 2 (K r ,-g) > . (75) 
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Figure 1: Left-hand: Time dependence of the expectation values ((-/Vi)t, (A^)*) for a field di- 
rection q = r = 1 and three different initial conditions. The dashed line shows the trajec- 
tory for (&i,&2) = ( 7r /2,0), the dashed-dotted line for {kx^k^) = (vr, 7r) and the solid line for 
(&i,&2) = ( 7r /2,— vr/2) . Right-hand: Trajectories ((a;), (y)) obtained by a direct integration of 
the time-dependent Schrodinger equation and the Hamiltonian (J^) with potential (|78[) for the 
same parameters and initial conditions as on the left-hand side. 

Here we also used that l-K^wl < 1 and < 2|iT Ujl) | 2 — l-K^u^l < 1 for a localised initial state. 
The second relation can be shown by using equation (|69|) and the relations for the ?9-function. 
Inserting this inequality into equation (|74jl shows that the quadratically increasing term \x[ r ' 9 ^ \ 2 
contributes to A 2 ^ (t) for every value of K u v resp. fcu,v Therefore, every Gaussian initial 
wavepacket shows a quadratic dispersion. 

The tight-binding results can be compared with a direct integration of the time-dependent 
Schrodinger equation for the Hamiltonian lfl|). The right-hand side of figure ^ shows the trajec- 
tories of the expectation values ((x)t, (y)t) f° r a Gaussian wavepacket 

il>(x,y) = (2^ 2 )- 1 /2 e -(x 2 +j/ 2 )/4 ( r 2 +ifc 1 a ; +ifc 2 y ( 76 ) 
with a = 107T which is projected onto the first Bloch-band: 

ipl(x, y) = (x, y\ipi) = ^2(m, n\rp) {x, y\m, n) (77) 

for three initial momenta (/ci,^) = (tt/2,0), (vr,7r), (vr/2, — vr/2). The system parameters are 
h = 2.81/27T and F = 0.1 in scaled units. The potential 

V(x,y) = cos(2-7T2;) + cos(2-7ry) — cos(27r:r) cos(2-7ry) (78) 

generates a lowest Bloch band in the energy dispersion relation with width Ai = A2 = 1/4 
along the axes and a width A3 = 0.008 along the main diagonals, so that the results can be 
directly compared with the tight-binding calculations on the left-hand side of figure ^ The small 
distortion compared to the tight-binding results can be explained by the couplings to further 
sites in the single-band model. 
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7 Concluding remarks 



In the present paper we have investigated Bloch oscillations in a two-dimensional tight-binding 
resp. single-band model. Using an algebraic approach, the evolution operator has been expressed 
in a product form for an arbitrary time dependent field. It has been shown analytically for the 
expectation values that an initial state in general does not perform a strictly periodic centre of 
mass Bloch oscillations in the single-band approximation. While the projection of this motion 
onto the field direction still oscillates, the projection orthogonal to the field additionally shows 
a directed motion. Only for special conditions that depend on the initial wavepacket and not 
on the system parameters the strict periodicity of the Bloch oscillations can be restored. The 
width of a wavepacket behaves similar: It shows a systematic quadratical dispersion orthogonal 
to the field that can only be suppressed for special initial wavepackets. 

Finally, the similarity to the one-dimensional tight-binding system with an additional har- 
monic driving, F(t) = Fq + F\ coswt should be emphasised. For an integer ratio between the 
Bloch frequency and the driving frequency, = (g G N) a wavepacket shows a systematic 
dispersion. But, in contrast to the two-dimensional dispersion , the dispersion can be suppressed 
by adjusting the field strength F\ and the driving frequency u> (dynamical localisation) [18-20]. 
Therefore, the two-dimensional tight-binding system with an additional time-periodic driving 
can be expected to show a multitude of interesting phenomena which have to be analysed in 
more detail. 
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